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We investigate possible phase transitions among the different quantum anomalous Hall (QAH) 
phases in a zigzag graphene ribbon under the influence of the exchange field. The effective tight- 
binding Hamiltonian for graphene is made up of the hopping term, the Kane-Mele and Rashba 
spin-orbit couplings as well as the Haldane orbital term. We find that the variation of the exchange 
field results in bulk gap-closing phenomena and phase transitions occur in the graphene system. If 
the Haldane orbital coupling is absent, the phase transition between the chiral (anti-chiral) edge 
state V — +2 [v = —2) and the pseudo-quantum spin Hall state {v = 0) takes place. Surprisingly, 
when the Haldane orbital coupling is taken into account, an intermediate QSH phase with two 
additional edge modes appears in between phases v = -\-2 and v = —2. This intermediate phase is 
therefore either the hyper-chiral edge state of high Chern number v — +4 or anti-hyper-chiral edge 
state oi V = —4 when the direction of exchange field is reversed. We present the band structures, 
edge state wave functions and current distributions of the different QAH phases in the system. We 
also report the critical exchange field values for the QAH phase transitions. 

PACS numbers: 71.70.Ej, 72.25.Dc, 73.43.Nq, 81.05.ue 



I. INTRODUCTION 

The anomalous integer quantum Hall effect observed 
in monolayer graphenes subjected to an external mag- 
netic field [l|, Q has recently attracted considerable at- 
tention. A theoretical investigation Q showed that 
plateaus located at the half-odd-integer position origi- 
nate from an additional Landau level at zero energy PjQj 
which is unlike the behavior of the conventional quan- 
tum Hall effect observed in two-dimensional (2D) het- 
erostructure semiconductors. @ The quantum Hall effect 
has also been experimentally observed in AB-stacked bi- 
layer graphenes 7], and this has been studied theoreti- 
cally as well. [8|, |9|] Furthermore, much experimental ev- 
idence is available for the existence of AA-stacked bi- 
layer graphenes. [l^ Interestingly, AA-stacked bilayer 
graphenes has been shown to exhibit zero transverse con- 
ductivity. [ll| 

According to Laughlin's gauge invariance argument, 
the sample edges are essential in generating the local- 
ized current-carrying states (edge states) [ij, • The 
edge states on the sample boundary are protected by the 
bulk band structure topology which is a manifestation 
of the Chern number, as elucidated by Thouless et al. 
(TKNN). fll The TKNN integer (or Chern number) 
relates the topological class of the bulk band structure to 
the number of chiral edge states on the sample boundary 
{bulk-boundary correspondence) and hence gives rise to 
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the quantized Hall conductivity axy = ve^/h. The pre- 
cise quantization of the Hall conductivity arises in the 2D 
electron system with an integer filling of the Landau lev- 
els. The Chern number corresponding to the number of 
the chiral edge currents equals to the number of Landau 
levels below the Fermi level. When the system undergoes 
a phase transition from one chiral edge state to another, 
the corresponding Chern number varies discontinuously 
from one integer u to v ±1. 

The Chern number must vanish in a system with time 
reversal symmetry (TRS). A TRS breaking mechanism is 
thus required for a 2D system to achieve a non-zero Chern 
number, either with or without the Landau levels. It has 
been shown that the chiral edge state in the quantum 
Hall phase is related to the parity anomaly of 2D Dirac 
fcrmions. [l^ [3 Therefore, in a remarkable paper [l7| . 
Haldane constructed a tight-binding Hamiltonian in the 
2D honeycomb lattice with a staggered magnetic field 
that produces zero average magnetic flux per unit cell 
(i.e., no Landau levels) and showed that the gapped state 
exhibits the quantum Hall phase with v = ±1. In this 
sense, the Haldane model is the prototype for the quan- 
tum anomalous Hall (QAH) effect. The relationship be- 
tween the Chern number and the winding number of the 
edge state was investigated in Ref. [l8| . 

On the other hand, when the bulk band gap of a sys- 
tem having a spin degree of freedom is opened due to the 
spin-orbit interaction, the system might be in the quan- 
tum spin Hall (QSH) state where the gapless edge states 
appea ring on the sample boundary are protected by the 
TRS. [19| The quantization of the spin Hall conductivity 
has been predicted in a graphene system with the Kane- 
Mele spin-orbit interaction as well as in a semiconductor 
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superlattice. [20l - |22j The quantization of the spin Hall 
conductivity, however, may be destroyed by the parity- 
breaking perturbations via spin non-conserving term or 
disorder. The associated topological invariant classifying 
the band structure topology of the time- reversal invariant 
systems is a Zi topological index [23,[2l|. The connection 
between the Chern number and the Zi topological index 
is explained in Ref. 23] . The Zi topological number rep- 
resents the number of the Kramer pairs of the gapless 
edge modes. An important result of this classification is 
that these gapless edge modes with an odd number of 
the Kramer pairs in the 2D systems [20l - [2^ and an odd 
number of surface Dirac cones in the three-dimensional 
(3D) systems [13, HH] are robust to impurity scattering; 
the other systems are just a conventional band insulator. 

Recently, another topological invariant (i.e., the spin 
Chern number) has been proposed by Sheng et al. [26|, 
and can be evaluated by imposing twisted boundary con- 
ditions on a finite sample. In Ref. [27| . it has been shown 
that the spin Chern number and Z2 topological orders 
would yield the same classification by investigating the 
bulk gap-closing phenomena in the time-reversal invari- 
ant systems. The phase diagram of the 3D QSH systems 
has been investigated systematically. [28| The topolog- 
ical winding number related to the spin edge states of 
graphene with the Kane-Mele Hamiltonian has also been 
studied. [29| Furthermore, the bulk-boundary correspon- 
dence is generalized to classify topological defects in insu- 
lators and superconductors, where the gapless boundary 
excitations are Majorana fermions. [30| 

In this paper, we first model the bulk graphene and 
also a zigzag graphene ribbon in the presence of the ex- 
change field [31] by using the Kane-Mele-Rashba Hamil- 
tonian [see Eq. ([!])]. Here the spin degeneracy is lifted 
by the TRS breaking term (i.e., the exchange field) and 
the z —z mirror symmetry is broken by the Rashba 
term. We calculate the Chern number of the bulk sys- 
tem as a function of the exchange field strength. Further- 
more, we study how the edge current in the corresponding 
graphene ribbon varies during a phase transition induced 
by the exchange field. 

The quantum anomalous Hall effect in Hgi_yMnyTe 
quantum wells (32| and tetradymite semiconductors 
(Bi2Te3, Bi2Se3, and Sb2Te3) [SS] has been investigated. 
Graphene with the Rashba spin-orbit coupling [a, see 
Eq. (Pel) ] and the exchange field has also been studied be- 
fore. 34] However, in Ref. Isij, the Kane-Mele spin-orbit 
couphng [A, see Eq. (pb)) ] was neglected because it was 
thought to be weaker than the Rashba spin-orbit interac- 
tion. In this paper, we find that, in the presence of both 
the Rashba and Kane-Mele couplings, a phase transition 
from either a chiral (i/ — +2) or anti-chiral (i/ = —2) edge 
state {i^ — ±2) to the pseudo-QSH state — 0) would 
occur in the graphene ribbon, because of the change of 
the Chern number due to the bulk gap-closing phenom- 
ena. This phase transition is different from the transition 
between the QSH state and the insulator state when the 
exchange field is absent. 



We then add the Haldane orbital coupling term which 
couple the electron orbital motion to the exchange 
field [13], to the Kane-Mele-Rashba Hamiltonian for 
graphene. Interestingly, we find that this leads to an 
anomalous change in the Chern number pattern. Note 
that the Haldane orbital term does not lift the spin de- 
generacy. Furthermore, we find that the presence of the 
Haldane orbital coupling would give rise to a new inter- 
mediate phase between phases ly = +2 and ly — —2. This 
intermediate phase has two new edge modes, and is thus 
either a hyper-chiral edge state with v — +4 or an anti- 
hyper chiral edge state with = — 4 when the direction 
of the exchange field is reversed. 

The rest of this paper is organized as follows. In SecHIl 
we describe the effective tight-binding Hamiltonian for 
graphene used in this work. In Sec. IIIIl we report the 
energy bands of a graphene ribbon in the presence of 
the exchange field. In Sec. IIVI we present the phase 
transition and the variation of the Chern number with 
the exchange field in the Kane-Mele-Rashba system. In 
particular, we show that the graphene ribbon undergoes 
a phase transition from the chiral (or anti-chiral) state 
to the pseudo-QSH state. In Sec. |Vl we show that a 
hyper-chiral (or anti-hyper-chiral) state would appear in 
between the chiral and anti-chiral states in the Haldane- 
Rashba system. The conclusions are given in Sec. IVIl 

II. EFFECTIVE TIGHT-BINDING 
HAMILTONIAN FOR GRAPHENE 

We consider the effective tight-binding model for 
graphene given by the Kane-Mele-Rashba Hamiltonian 

Hkmr — Ht + H\ + Ha , (1) 

with 

Ht = tJ2 4ci^ (2a) 
<»,j> 

Hx = tX ^ 4^.i^:jC,, (2b) 

Ha^ia ^ 4(s X dij)^Cj. (2c) 
<i-j> 

The symbols < i,j > and <C i,j ^ denote the nearest 
neighbors and the next nearest neighbors, respectively. 
The Hamiltonian Ht is the tight-binding energy for the 
nearest-neighbor hopping. The Kane-Mele Hamiltonian 
Hx describes the intrinsic spin-orbit interaction. The 
site-dependent Haldane phase factor ^ Vij is defined 
as 



where denotes the vector from one carbon atom to 
one of its nearest neighbors. Two vectors di and d2 are 
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required to represent the second neighbor hopping (see 
Fig.[T]). In the two-dimensional case, the non-zero compo- 
nent becomes a sign function and we take the values 
of ±1 (i.e., counterclockwise/clockwise). The extrinsic 
spin-orbit interaction is described by the Rashba Hamil- 
tonian Ha, which can be produced by, e.g., applying an 
electric field E perpendicular to the graphene sheet. Ha 
is proportional to E • (s x d^j), where E = E'^Cz and dij 
denotes the vector from site i to site j (see Fig. [1]). 

Recent ah initio density functional calculations showed 
that intrinsic ferromagnetism in pure and on-top-Fe- 
doped graphene monolayers may exist, [s^ . [35| Further- 
more, proximity-induced ferromagnetism in graphene 
was recently reported. [36| Therefore, we consider the 
interaction of the 2D electrons in graphene with the ex- 
change field produced by the ferromagnetism. [3l| The 
coupling of the orbital motion and also spin of the elec- 
trons on graphene to the exchange field would give rise 
to an additional Hamiltonian: 



with 



Hex — Hj + Hfj , 
Hj = 7^4szQ, 

i 



(4) 

(5a) 
(5b) 



where 7 is the (rescaled) exchange field strength. The 
coupling 7 is proportional to Jcff/i^, where Jeff is the ex- 
change interaction and /^'j, is the effective magnetic mo- 
ment associated with the exchange field. The magnetic 
field generated by ji'^ is denoted as 7'. The Hamiltonian 
iJ-y describes the response of an electron spin magnetic 
moment to the exchange field a la Zeeman effect. 

In the meantime, the orbital angular momentum of an 
electron in graphene would be coupled to the exchange 
field because of its associated orbital magnetic moment. 
The Haldane phase factor Vij behaves like an effective 
orbital angular momentum, and hence gives rise to the 
interaction between the electron orbital motion and the 
magnetic field 7', as described by Eq. ((5b| . where is 
a function of 7'. Spatial parity symmetry requires /3(7') 
to be an odd function of 7'. Unlike the interaction be- 
tween the spin and exchange field, the energy of the Hal- 
dane orbital coupling cannot be linear in the exchange 
filed 7'. Instead, the response of Haldane orbital motion 
would be saturated rapidly because the exchange field 7' 
alters the orbital velocity of electrons and induces an or- 
bital magnetic moment against it. Phenomenologically, 
we can adopt the simple yet sensible approximation: 



/3(7')«/?sgn(7), 



(6) 



where we use 7 instead of 7' for simplicity since the sign 
function is independent of the field strength but its direc- 
tion. In the present study, the sign of Jeff is fixed, and 
hence the sign change of 7 corresponds to the change 




FIG. 1: (Color online) (a) A segment of a zigzag graphene 
ribbon with its unit cell marked by the red dashed lines, (b) 
Illustrations of Haldane phase factors v^j, vectors Aij, bulk 
basis vectors (01,02) and bond length a. (c) The first Bril- 
louin zones of bulk graphene (gray region) and the zigzag 
graphene ribbon (red double arrow) in the 2D fc-space. 



in the direction of /i^, which is experimentally possible. 
Accordingly, we choose the constant /3 to be negative to 
have a diamagnetic response to the magnetic field 7'. 



III. CHERN NUMBERS AND EDGE CURRENT 
CHIRALITY 

The total Hamiltonian for graphene in the presence of 
the exchange field is given by i? = i?KMR + -ffcx- For 
the bulk graphene, the Hamiltonian i?(k) which satisfies 
the periodicity i?(k) = i?(k + G) (G stands for a 2D 
reciprocal-lattice vector) , is given by 



iJ(k) = 



/ A'Z + 7 

X~iY 


V ~iaM*_ 



X-l-iY 
-A'Z-K7 
iaM+ 






-iaM\ 

- 7 
X-iY 





X + iY 
A"Z - 7 



(7) 

where A' = A + /3sgn(7) and A" = A — /3sgn(7). The state 

vector is represented by tl)"^ = (c^^i, Cks-|-, Cksj.), 
where A and B denote the two different sublattice points 
in the unit cell, respectively, and the arrows represent 
the spin directions. The matrix elements are given by 
X = t[l -I- 2cos(fc^)cos(3fc;)], Y = t[2cos(fci)sin(3fc;)], 
Z = 2sin(2fc;)-4sin(fc;)cos(3fc(,), M+ = [-l+2cos(fc;- 
f ) cos(3/c' )] + i[2 cos(fc; - f ) sin(3fc' )], and M_ = [-1 
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2cos(fc; + f )cos(3fc^)] +i[2cos(fc^ + f)sin(3fcp], where 
the two variables k'^ and k'y are defined as k'^ = ^k^a 
and k'„ = %a, respectively. Note that along the ky — 



profile, the two points fc^ = ±^ are just the K and 
if' points in the Brillouin zone of bulk graphene [see 
Fig. [UJc)], respectively. After the eigenvalue equation 
iJ(k)IV'nk) = i^raklV'rak) IS solved, the Berry curvature 
{^^xy ) of the nth band can be calculated using 



E 



2 lm(-0„k I V'n'k) (V'n' k I «y I V'ti 



The Chern number is then obtained by summing the 
Berry curvatures VLxy for all the occupied states below 
the Fermi level for each k and subsequently integrating 
over the entire first Brillouin zone: 



dk,dkyni-\k). 

„ Jbz 



(9) 



The bulk Hamiltonian Eq. ([7]) is simplified greatly if 
we consider the following simple systems: 



1. Kane-Mele system: H] 



KM 



Ht + H\ + 



2. Rashba system: H-r — Ht + Ha + H^ 



3. Haldane system: H^ = Ht + H, 



+ H^. 



The sign of Chern number indicates the chirality of the 
edge current. To verify the occurrence of the edge cur- 
rents, we compute the energy band structure for a zigzag 
graphene ribbon. The unit cell of the zigzag graphene 
ribbon is shown in Fig. [Ija), where the ribbon direction 
is denoted by the x axis and the transverse direction is 
along the y direction. The width of the zigzag ribbon 
(W) is 75 a, where a is the bond length [see Fig. [I{b)], 

1. e., there are + 1 = 101 C atoms in the transverse 
direction [see Fig. [IJa)]. The nearest neighbor hopping 
integral t = 1. 

Figure 2 shows the ribbon band structure and the edge 
state probability distribution in the Kane-Mele system 
(a), the Rashba system (b), and the Haldane system (c). 
The Fermi level is assumed to be above zero, as indicated 
by the dashed horizontal line, and thus, has four inter- 
sections with the conduction bands, denoted as A, B, C, 
and D, in the left panels in Fig. 2. This gives rise to 
four edge currents on the ribbon edges, as indicated by 
the A, B, C, and D arrows in the right panels in Fig. 

2. The direction of an edge current, denoted by an ar- 
row, is given by / = — [elu^: where the electron group 
velocity is determined using — dE\^/dkx. The A and 
B states have the same velocity direction, which is op- 
posite to that of the C and D states. Hereafter, we use 
the notation (J^, Ir) to express the charge current distri- 
butions on the left-hand side and right-hand side edges, 
respectively. In terms of the bulk-boundary correspon- 
dence, for each of the three systems, the pair A and D 
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FIG. 2: (Color online) Calculated energy bands (left pan- 
els), edge state probability (middle panels) and charge cur- 
rent (right panels) distributions in the zigzag ribbon in the 
presence of the exchange field. The Fermi level (the dashed 
line in the left panels) Ef ~ 0.05t. (a) Kane-Mele system 
(A = 0.06t and 7 = 0.2t). (b) Rashba system (a = 0.2t and 
7 = 0.2t). (c) Haldane system (/3 = -0.07t and 7 = 0.2t). 



would form a single handed loop (the turning point is at 
infinity in the a; direction), and the pair B and C would 
constitute the other loop of the opposite handedness, as 
can be seen from the probability distribution shown in 
the middle panels in Fig. [2] 

In the Kane-Mele system Hkm, the current distribu- 
tion is {Ibd, Iac), as indicated in the right panel in 
Fig. dja). The two edge states A and C are on the 
same edge, and so are the B and D states. As men- 
tioned above, the handedness of the current loop due to 
the A and D edge states would produce a Chern num- 
ber of —1 while that of the pair B and C would give 
a Chern number of -1-1. Therefore, the Kane-Mele sys- 
tem is composed of two integer quantum Hall subsys- 
tems, namely, (i^ — +1) ffi (j/ = —1), [13, 113 and has 
1/ = (-1-1) -I- (—1) = 0. Since this state has the same 
distribution of the edge currents as that of the quantum 
spin Hall case with the TRS [l^l, except that the TRS 
is broken here, we call this state as the pseudo-quantum 
spin Hall state. 

In the Rashba system iJpj , the current distribution is 
{Iab, Icd), as shown in Fig. [DJb), which constitute a 
paramagnetic response to the exchange field. Both I a 
and Ib are located at the same edge, confirming that the 
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Rashba system has a Chern number of +2, since the two 
edge current pairs have the same chirahty. 3J] It is im- 
portant to note that both the Chern number and the cur- 
rent distribution {Iab, Icd) in the Rashba system is in- 
variant under the transformation a — > —a. On the other 
hand, the current distribution becomes {Icd, Iab) when 
the direction of the exchange field is reversed. Therefore, 
the Rashba system is equivalent to two integer quantum 
Hall subsystems, namely, (j/ = -|-l)©(i/ = +l)for7>0 
or {f = -1) ® (i/ = -1) for 7 < 0. m 

In the Haldane system TJh, the current distribution is 
{Icd, Iab), as shown in Fig. [2ljc). Both I a and Ib are 
also located at the same edge, but the chirality of the 
edge current is opposite to that of the Rashba system, as 
a result that the Haldane system with /? < exhibits a 
diamagnetic response to 7'. The Chern number of this 
system is = —2. Therefore, the Haldane system, being 
diamagnetic, is equivalent to two integer quantum Hall 
subsystems, namely, (v = —1) (B = —1) for 7 > or 
{i' = +1) ffi (i/ = +1) for 7 < 0. 

In the next two sections, we will consider the following 
two combinations of the three simple systems discussed 
in this section: 



1. Kane-Mele-Rashba system: Hi = Hkmr 

2. Haldane- Rashba system: II2 = Ht +Hp + Ha 



We find that, because of the bulk gap-closing phenomena, 
both systems will undergo a change of the edge current 
chirality caused by varying the exchange field. 



IV. PHASE TRANSITION IN THE 
KANE-MELE-RASHBA SYSTEM 



In this section, we will neglect the Haldane orbital cou- 
pling term of Eq. (jSbl) . We will find that the phase tran- 
sition is different from the QSH phase transition in the 
presence of the exchange field. We consider the interplay 



between and iJi 



KMR- 



Hi — Hkmr + 

^ Ht + Ha + H\ + H~, 



(10) 



In the presence of both Kane-Mele and Rashba spin-orbit 
couplings, the phase transition between the chiral (or 
anti-chiral) state and the pseudo-QSH state must occur 
when the bulk gap-closing phenomena take place. On 
the other hand, although the locations of four currents 
{Ibd, Iac) become {Iac, Ibd) under the transformation 
A — > —A, the phase transition between (anti-) chiral and 
pseudo-QSH states still applies because both Ibd and 
Iac correspond to 1^ = 0. 

In order to verify the the phase transition in the finite 
system (the graphene ribbon), we use the expectation 
value of position y (i.e., (y)) as a parameter for specifying 
the angular momentum of the current in the system, and 
define (y) = {y}A + {y}B- When the Kane-Mele couphng 
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FIG. 3: (Color online) Expectation value (y) and Chern num- 
ber as a function of "//t. (a) The Kane-Mele-Rashba system 
with A = 0.06t, a = 0.05t, and 7 ranging from — 0.5t to 0.5t. 
(b) The Haldane- Rashba system with /3 = — O.lt, a — 0.5t 
and 7 ranging from —t to t. 



A is dominant (I7I < 7*^, see below), {y)A and {y)B are 
on the opposite sides of the ribbon, and thus (y) /W = 0. 
When the Rashba coupling is dominant (I7I > -f'^, see 
below), {y)A and (y) b are on the same side of the ribbon. 
The quantity (y) /W in the Rashba dominant system, 
however, would not reach a saturated value {y)/W — ±1 
owing to the finite size effect as the perfect edge states 
is obtained only when the ribbon width W is infinite. 
During the phase transition, the wave function starts to 
mix with each other in the central region of the ribbon, 
and thus, (y) is expected to deviate from either or ±1. 

Let us consider the case of A = 0.06i, a = O.Q5t and 
7 ranging from — 0.5i to 0.5i in Eq. ([T0|) . as an example. 
For 7 > 0, we find that (y) decreases to zero near some 
magnitude of 7 (see below) and does not change sign, as 
shown in Fig. [Sja) . The pattern of (y) in the 7 < region 
is the parity symmetry of that in 7 > 0. We find that the 
expectation value (y) changes sign when the direction of 
the chiral current is reversed. 

Based on the bulk-boundary correspondence [13 > 
the existence of the phase transition is supported by eval- 
uating the critical values of the exchange field. The crit- 
ical value of the exchange field (7^^) for the occurrence 
of the phase transition is determined by the bulk gap- 
closing phenomena, where the bottom of the bulk con- 
duction band (Ec) and the top of the bulk valence band 
(Ey) become degenerate, namely, Ec — Ey = at 
It can be shown that the degenerate point is located at 
k'^ = The critical value for the exchange field is 
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^fiven by 



± 



4A 



(11) 



which is obtained for a non-zero Kane-Mele couphng that 
satisfies 7/A < 3 a/3 and a/A < v^- The presence of A 
causes the critical value for the exchange field to shift 
from 7"^ = (Rashba system) to a non-zero value. The 
magnitude 7*^ corresponds to the location where the bulk 
valence and conduction bands become degenerate and 
the Chern number starts to jump from one integer to 
another. 

For a system with given a and A, if 7 > 7^ (7 < 7^), 
the system is in the chiral current state (pseudo-QSH 
state). When 7 = -l-|7^| (or 7 = — 17'^|), the conduction 
and valance bands touch at k'^ — 27r/3 and k'^ = — 27r/3 
simultaneously. Because the I a current direction is op- 
posite to Id , the exchange of the locations of I a and Id 
results in a change of chirality [see Fig. [Hfa)]. The cor- 
responding variation of the Chern number is also shown 
in Fig. |3^a). Therefore, we find that the Chern number 
jumps from 1^ = to v — ±2. In this case, the critical 
values 7'^ — ±0.293i are in agreement with the numerical 
result. 



V. PHASE TRANSITION IN THE 
HALDANE-RASHBA SYSTEM 

In this section, we will neglect the Kane-Mele coupling 
A. We will show that the presence of the Hamiltonian 
iJ^ creates two new edge modes between the two 1/ = ±2 
phases. Interestingly, these intermediate states are either 
the hyper-chiral state {v = +4) or anti-hyper-chiral state 
(1/ = —4). The Hamiltonian is given by 



II2 "111+ Ha + Hp + H^ 



(12) 



As described in Sec. IIIIl the Haldane system has v = —2 
and the Rashba system has v — +2 when the exchange 
field is positive. If the phase transition occurs in this sys- 
tem, the bulk gap-closing phenomena must take place. In 
the following, we show that the Haldane-Rashba system 
has two critical values of the exchange field. 

For the sake of discussion, we consider the region with 
7 > 0, and focus on the region A:^ > because the be- 
havior of the corresponding degenerate points in fc^ < 
is the mirror symmetry of that in fc^ > 0. In the pres- 
ence of a only, there are two degenerate points. One is 
at k'^i = 2tt/3 (i.e., K point), and the other is at 



k' 



/ 2a^-l 
' I 2a2 + 2 



(13) 



Eq. p3)) shows that the degenerate point depends on the 
strength of the coupling and thus varies with the mag- 
nitude of a. However, the two degenerate points appear 
simultaneously, namely, there is only one critical value 



for the exchange field, 7"^ = 0, even when the coupling A 
is considered [as shown in Fig.[31[a)]. 

Unlike the Kane-Mele coupling A, we find that the cou- 
pling l3 will close the bulk gap twice at two different mag- 
nitudes of the exchange field. Namely, there are two de- 
generate points appearing at two different magnitudes of 
the exchange field. Therefore, the Chern number varies 
discontinuously from v = +2 to v = —2 through an in- 
termediate state. One of the degenerate points is at the 
K point; the corresponding critical value of the exchange 
field is 



-SVS/S ; kl\ = 27r/3. 



(14) 



However, the second degenerate point for the Hamilto- 
nian Eq. is different from k'^2 expressed in Eq. (|13l) . 
The second degenerate point is determined by the condi- 
tion Ey — Ec = and can be expressed as 
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-4V3/3(1-F2) ; fc^'2 = cos-i(F), (15) 



where F satisfies the following equation: 

-1 + 2^2+32/32 (1 - + F^) -2F {l + a^ + 16/3^) = 0. 

(16) 

Numerically, Eq. can be solved for a given set of a 
and l3. When /3 = 0, Eq. ^ is the solution of Eq. (dH) 
and the critical value of the exchange field is 7f = 7f = 
0, which is in agreement with that in the case of the 
Rashba system. For the present case, a — 0.5t and /? = 
—O.lt. The two critical points are 7^ — 0.5196i and 72 = 
0.68995t. Therefore, the conduction and valence bands 
first touch at 7 = 7^, and the bulk gap would re-open 
when 7i < 7 < 7f (referred to as the intermediate state) . 
The two bands would touch the second time at 7 = 72. 
The bulk gap is open when 7 > 72. The calculated Chern 
number as a function of ^/t is shown in Fig.[21[b). 

Surprisingly, the intermediate state between the crit- 
ical values 7f and 72 shows v = —4. The band struc- 
ture and current distribution of the intermediate state 
are shown in Fig. We find that in the presence of 
Haldane orbital effect, the system establishes two new 
edge modes: one is the pair A2 and D2, and the other 
is the pair C'2 and i?2- Furthermore, the current distri- 
bution II (and I^) also shows that the four currents in 
II have the same chirality. Importantly, we find that un- 
like the quantum Hall plateau, the Chern number is not 
restricted in changing from one integer v to the next in- 
teger 1. Instead, a higher Chern number can exist in 
a system with a spin-orbit interaction if the orbital effect 
is also taken into account. 

The Chern number [see Eq. (j9])] can be written as 
u = dk'^dkyflxy. The bulk band structure and the 
corresponding Berry curvature flxy along the ky = pro- 
file are shown in Fig. |3Jc) and (d) , respectively, where we 
use a — 0.5t, fi = O.li and 7 = 0.6t. It can be shown 
that if the linear term /3i7' is taken into account, the 
Chern number is still —4. This clearly shows that the 
anomalous Chern number v = —4 is due to the second 



7 




-W/2 y W/2 -W/2 y W/2 



FIG. 4: (Color online) The band structure (a), the edge state 
probability and current distributions (b), bulk energy bands 
(c) and Berry curvature (d) of the intermediate state in 7^ < 
7 < 72 (a = 0-5t, 13 = -Q.lt, and 7 = 0.6t) (see text) in 
the graphene ribbon. The Fermi level Ef is at O.Olt [the red 
dashed line in (a)]. In (a), the 8 edge states are marked as 
A, B,C , D, A2, B2,C2, and D2, respectively. In (c) and (d), 
the bulk energy bands and Berry curvature are plotted along 
the ky = profile. 

nearest-neighbor hopping in graphene. We believe that if 
the third nearest-neighbor hopping is considered in Hal- 
dance orbital effect and spin-orbit interaction, a higher 
Chern integer of, e.g., —6, may be obtained. 

Let us define (y) = {y)At + (y)Bi + (j/)a2 + {y)B2 for 
the intermediate state. The calculated variation of {y) 
with 7 is shown in Fig. [2Ib). Apart from the occur- 
rence of the intermediate state, the expectation value 
{y) = {y)A + {y)B changes sign as the exchange field 
is swept through the phase transition, and this is accom- 
panied by a change of the Chern number from ly = +2 to 
v = —2, as shown in Fig. Eljb). The Chern number ob- 
tained from the bulk Hamiltonian represents the number 
of the perfect edge states. Note, however, that because 
of the finite size effect, {y)/W cannot reach the satu- 
rated value {y)/W = ±2 [see Fig. [Sjb)]. Interestingly, 
we find that it is not necessary to reverse the direction 
of the exchange field in order to flip the current chirality 



in this case. Therefore, the graphene can be brought to 
either the paramagnetic phase or the diamagnetic phase 
by adjusting the magnitude of the exchange field. 

Very recently, Tse et al. [Ill also proposed that the 
Hall conductance can be quantized as Uxy = 4e^//i, 
albeit^ in bilayer graphene with the Rahsba coupling 
under the influence of an external gate voltage. In 
the present work, in contrast, we show that in single- 
layer grapnene, the quantized Hall conductance can 
be axy = 4e^/ft. when the Haldane orbital effect is 
considered. Furthermore, we find that the change in the 
quantized Hall conductance can be achieved by varying 
the exchange field instead. 



VI. CONCLUSIONS 

In summary, we find that the edge current chirality in 
a graphene ribbon can be flipped by varying the exchange 
field. The resultant phase transition of the current chi- 
rality is caused by the bulk gap-closing phenomena; that 
is, the phase transition is due to the topological effect of 
the bulk band structure of graphene. We show that the 
paramagnetic response in the Rashba system can exhibit 
ether the chiral or anti-chiral state, and thus, the Hall 
conductance is quantized as axy — ±2e^/ft.. Wc find that 
the Kane-Mele system has the pseudo-QSH state. For 
the Kane-Mele-Rashba system, the transition between 
the chiral (or anti-chiral) current state and the pseudo- 
QSH state can be achieved by varying the strength of the 
exchange field. 

Unlike the Rashba system, the Haldane system ex- 
hibits a diamagnetic response to the exchange field, and 
the quantized Hall conductivity is axy = ±2e^/ft.. How- 
ever, the competition between a and /? leads to a phase 
transition between the diamagnetic and paramagnetic re- 
sponses, and hence an intermediate phase. Interestingly, 
this intermediate phase has two new edge modes and 
is thus a new quantum anomalous Hall state with high 
Chern number u = ±4. The corresponding quantized 
Hall conductance is axy = ±4e^//i in the graphene rib- 
bon in the absence of Landau levels. 
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